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Recent quantum Hall experiments conducted on disordered graphene pn junction provide evidence 
that the junction resistance could be described by a simple Ohmic sum of the n and p mediums' 
resistances. However in the ballistic limit, theory predicts the existence of chirality-dependent 
quantum Hall plateaus in a pn junction. We show that two distinctively separate processes are 
required for this ballistic-Ohmic plateau transition, namely (i) hole/electron Landau states mixing 
and (ii) valley isospin dilution of the incident Landau edge state. These conclusions are obtained 
by a simple scattering theory argument, and confirmed numerically by performing ensembles of 
quantum magneto-transport calculations on a 0.1/xra-wide disordered graphene pn junction within 
the tight-binding model. The former process is achieved by pn interface roughness, where a pn 
interface disorder with a root-mean-square roughness of lOnra was found to suffice under typical 
experimental conditions. The latter process is mediated by extrinsic edge roughness for an armchair 
edge ribbon and by intrinsic localized intervalley scattering centers at the edge of the pn interface 
for a zigzag ribbon. In light of these results, we also examine why higher Ohmic type plateaus are 
less likely to be observable in experiments. 

PACS numbers: 



I. INTRODUCTION 

Graphene is a two-dimensional sheet of carbon atoms 
arranged in a honeycomb lattice with unique electronic 
properties i.e. linear energy dispersion with zero bandgap 
described by the relativistic Dirac equation [TJ [2] . This 
attribute manifests itself as anomaly in the quantum 
Hall regime [3j [4], where the Landau filling factor, 
goes by the non-conventional sequence of 4n + 2 and 
n = 0, 1, 2, . . . labels the Landau levels (LL). This leads 
to conductance plateaus at a = z/, expressed in units of 
the von Klitzing constant e 2 /h. Room temperature ob- 
servations [4] of these quantum Hall plateaus make it a 
potential state variable for possible future device applica- 
tions, such as quantum computation using Landau levels 
as the basic qubit [5 and pure spin current switching 
devices [6]. 

Recently, experiments by Williams et al. [7 on the 
magneto-transport of graphene pn junction in the quan- 
tum Hall regime found that the two terminal junction 
conductance (a pn ) exhibits new plateaus. These conduc- 
tance plateaus were predicted to follow a simple Ohmic 
conductance rule, first proposed by Abanin et al. [8], 
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where v n / p are the Landau levels filling factors in the n/p 
regions. The last expression in Eq. [T] presents the physics 
more lucidly, with (. . .) eq embodying the Landau modes 
mixing process along the pn interface for the electron 
current. Experiments [7], however, did not reproduce all 
the predicted plateaus, especially for the higher filling 



factor combination values i.e. such as (v n ^ p ) = (2,6) 
and (6, 6) [54 . This observation leads us in formulating 
our first question, "what degree of interface disorder is 
required to observe complete modes mixing, especially 
for the cases with higher combination values?" . 

Interestingly, the appearance of quantum Hall plateaus 
in a graphene pn junction is not isolated to only disor- 
dered samples. Tworzydlo and co-workers [9] showed the- 
oretically that new plateaus should be observed in per- 
fectly clean graphene ribbons with perfect edges. These 
plateaus are independent of (y n ,v p \ and depend only on 
the ribbons' chirality. We herein denote these ballistic 
plateaus as <j pn , and they are given by [9], 
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where the widths dictate the above possible outcomes. 
The underlying physics for this width dependence is ex- 
plained in terms of the valley isospin (for armchair-edge 
[9]) and parity (for zigzag-edge [10]) of the lowest lying 
Landau modes. The above revelation begs the second 
question, "Is pn interface disorder alone sufficient in in- 
ducing the plateau transition from a pn to a pn , or is edge 
disorder (or other intervalley scattering processes) also 
necessary?" . 

Theoretical studies on the effect of disorder on the 
magneto-transport properties of a graphene pn junction 
are few [9j [TTJ [12] . In the classical integer quantum Hall 
problem, the tight-binding Hamiltonian model [13] and 
the Chalker/Coddington's network model [14] are popu- 
lar approaches for studying magneto-transport in disor- 
dered system [15] . For this problem the former is more 
suitable, since the formalism inherently captures both 
quantum mechanical and atomistic effects which are be- 
lieved to play an important role. Essentially, one seeks 
the direct solution to the one-electron Schrodinger equa- 
tion, where the open boundary scattering problem [16] 
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is usually conceptualized within the Landauer-Biittiker 
quantum transmission point-of-view [TTl [18] . The quan- 
tum transmission function can be calculated using the 
Green's function [T9l [20] or wave function approach [21], 
the former being a more popular technique in recent 
years. Within this theoretical framework, Long and co- 
workers [11] conducted an intensive statistical study of 
disordered graphene pn junctions in the quantum Hall 
regime. Graphene ribbons with zigzag edges are consid- 
ered in their work. Bulk disorder was incorporated into 
their Hamiltonian in the form of on-site energy fluctua- 
tions. They found that the Ohmic type quantum Hall 
plateaus, <r pn , emerges with sufficient disorder strength. 
Similar conclusions had also been reached by Li and Shen 
[12], where pn interface disorder was considered for a 
zigzag ribbon, also in the form of on-site energy fluctu- 
ations. Both these numerical studies [TTJ [12] employed 
short range disorder potential which varies on the scale of 
the lattice constant. Although the actual sources of dis- 
order varies across different experimental samples [22] . 
disorder such as pn interface and edge roughness are not 
justifiably captured by short range disorder potential. 
Moreover, short range disorder potential serve to masked 
the effect of valley isospin, which is inevitably washed out 
by such disorder [12 . For example, Tworzydlo et al. [9] 
had employed a long range disorder potential, and found 
that the junction magneto-conductance of a zigzag rib- 
bon is extremely sensitive to the disorder, while an arm- 
chair ribbon shows the opposite behavior. 

The purpose of this paper is to quantitatively model 
the effect of pn interface and edge disorders and to estab- 
lish a coherent, conceptual framework for understanding 
the plateau evolution from & pn to a pn . We reason that 
one can distinguish two underlying mechanisms which are 
prerequisites for this ballistic-Ohmic plateau transition, 
namely (i) hole and electron Landau states mixing and 
(ii) valley isospin dilution. The former can be achieved 
via pn interface roughness, a long range potential type 
disorder. The latter process is mediated by extrinsic edge 
roughness for an armchair edge ribbon and by intrinsic lo- 
calized intervalley scattering centers at the edge of the pn 
interface for a zigzag ribbon. By performing ensembles 
of quantum magneto-transport calculations on a 0.1/im- 
wide graphene pn junction with pn interface (ID) and 
edge disorder (ED), we illustrate how the plateaus evolve 
from a pn to a pn for both armchair and zigzag ribbons. 

This paper is organized as follows. Section II be- 
gins with an introduction of our quantum transport 
model based on the Green's function and tight-binding 
approaches. Computational and numerical aspects are 
highlighted in this section, including the statistical mod- 
eling of the interface and edge disorder. Section III dis- 
cusses the magneto-transport across a pn junction in the 
ballistic limit. The underlying physics for the existence of 
the ballistic quantum Hall plateaus a pn are reviewed and 
discussed. Section IV studies the ballistic-Ohmic quan- 
tum Hall plateau transition in the presence of ID/ED 
for both armchair and zigzag ribbons. An accompany- 
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FIG. 1: (a) Illustrations of armchair, zigzag and anti-zigzag edge 
ribbons. The carbon layer numbering convention for an armchair 
edge ribbon used in this work is also depicted, (b) Schematic of 
the two terminal pn junction simulated in this work, illustrating 
also the ID and ED. 



ing simple Chalker/Coddington type scattering theory is 
presented to elucidate the underlying physics. Section V 
pertains to the analysis and comparison with recent ex- 
perimental data in the literatures followed by a summary 
of this work. 



II. QUANTUM TRANSPORT MODEL 

In quantum Hall experiments, a magnetic field of the 
order of 10T is common. This corresponds to a magnetic 
length of £b ~ 20nm, which approximates the spatial ex- 
tent of the ground state LL wavefunction. The device's 
width has to be appropriately chosen so that the LL's 
wavefunctions along the opposite edges do not overlap, 
especially for the higher LLs. However, finite computa- 
tional resources set a practical limit on the matrix size of 
the system Hamiltonian. Based on the above consider- 
ations, we employed a lOOnm-wide graphene ribbon pn 
junction device, which is schematically shown in FigjTJ^. 
Ribbons of both armchair and zigzag- type edges are con- 
sidered in this work [22] [23] , as shown in FigjTJi. We seek 
to compute the junction conductance (cr pn ) in the pres- 
ence of pn interface and edge disorder. For pedagogical 
purposes, we shall briefly review the quantum transport 
model employed in this work. 



A. Green's function theory with a tight-binding 
model Hamiltonian 



In this work, the device is described by tight-binding 
Hamiltonian given by 

H = Vid\ai + \Uj | exp (ij^l A • d\ j a\aj , (3) 

i ij \ t / 

where a\/di are the creation/destruction operators at 
each atomic site z, and V{ and Uj are the on-site poten- 
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tial energy and hopping energies [24 , 25 j. The simple one 
^-orbital description as described by Eq. [3] is sufficient 
for the modeling of the relevant energy bands of interest 
for electronic transport properties [55] . We assumed that 
\Uj \ = 3eV. In the presence of a perpendicular magnetic 
field B, tij will incorporate a Peierls phase <pij (A), where 
A is the corresponding vector potential. The phase (j)^ 
is assigned such that the magnetic flux through an arbi- 
trary area S satisfies the following Stokes law [26] , 



(4) 



where Q is the boundary of S. The parameter V{ is de- 
termined by the top/bottom gate electrostatics, which is 
known a priori. The effect of pn interface and edge dis- 
order are also described through V{ and respectively 
(discussions defered to next section). 

The Landauer-Biittiker approach pictures a device in 
which dissipative processes are absent but coupled to per- 
fect thermodynamic systems known as "reservoirs" . This 
approach has been very successful in modeling physical 
effects in a myraid of problems in the field of mesoscopic 
physics p~9j [20] [27] [28] . In numerical implementation, H 
is divided into "device" (Hd) and "contacts" regions. H d 
is constructed so that it captures the scattering region of 
interests, hence it is of finite matrix size. The contacts 
represent the semi-infinite regions of H, which character- 
ized the open boundary nature of the transport problem. 
Through simple algebras as detailed in [19], one could 
write a Green's function for the "device" as follows, 



G(e f ) = (e f -H d -Z s -X d y 



(5) 



where €f is the Fermi energy. £ s /d are conveniently 
known as the contact self-energies (subscript s/d for 
source/drain respectively), which could be expressed as 
^s/d = r s / d g a / d rl /d , where r s/d describes the coupling 
of the s/d contacts to the device and g s j d is the surface 
Green's function of the respective contacts. In this work, 
g s / d is obtained using an efficient iterative scheme out- 
lined in [29 . Direct matrix inversion of Eq. |5| proves 
to be computationally prohibitive. Therefore, one com- 
monly resorts to recursive type techniques, such as the re- 
cursive Green's function approach [30j [31], the renormal- 
ization method [32], or combination of both techniques 
[33] . After solving for G(e/), we can compute the device 
conductance at ef via (in units of ^-) [19], 



a pn = 2Tr (r s GT d G^) 



(6) 



where T s / d are known as the contact broadening func- 
tions which can be obtained from the respective self- 
energy i.e. T s/d = i(Yl s /d - S^ /d ). One can view Eq. 
[6] as just a different form of the Fisher-Lee expression 
[34 [19 . Other physical observable quantities such as the 
local density of states (LDOS), charge density (n(x,y)) 
and current density (j(x,y)) can also be obtained [33] . 
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FIG. 2: Intensity plot of the longitudinal current component for 
magneto-transport across an armchair graphene pn junction at 
B = 10T. We consider ground state Landau level injected from 
the left. The biased condition is such that the filling factor com- 
binations is (i/ n , Up) = (2, 6), where spin degeneracy accounted for. 
The depletion width is assumed to be 25nm and perfect edges is 
assumed. We plotted for the case (a) without ID and (b) with ID 
for ribbon with 401 carbon layers along the width. Similary, for 
ribbon with 400 carbon layers in (c) and (d). 



B. Modeling interface/edge disorders 

In this work, we are interested in two types of disor- 
der - interface and edge disorder (ID/ED). It is generally 
accepted that ID plays an important role in facilitating 
the mixing of the electron/hole Landau modes along the 
pn interface [8]. We define the pn interface as the equi- 
energy line of Dirac points separating the n/p regions. 
This interface will be highly susceptible to electrostatic 
influences of impurities due to ineffective charge screen- 
ing in the depletion region. The presence of electron- 
hole puddles due to trapped impurities in the oxide lay- 
ers [35] (or ripples [36]) was recently observed by Mar- 
tin and co-workers [37]. Charge density fluctuations of 
±1 x 10 11 cm -2 was reported in their work. This implies 
significant ID. Edge roughness was also recently charac- 
terized by Gupta et al. [38], where a root-mean-square 



4 



roughness of about 3nm was observed in micromechani- 
cally produced ribbons. 

In order to model ID and ED, we devise a simple algo- 
rithm for generating a one dimensional roughness profile, 
I(y). I(y) can be expressed in its Fourier components, 

I(y) = Y J A n sin{^) (7) 

n 

where W is the device width, and A n is the amplitude 
for the sine components given by, 

A n = R(D 1 )exp(-%) (8) 

where R(Di) outputs a uniformly distributed random 
number between ±Di. Eq. [8] represents the power spec- 
trum of the roughness morphology, however there is cur- 
rently no characterization of ID/ED morphologies to jus- 
tify such assumption. Nevertheless, it is known that an 



exponential power spectrum describes the surface mor- 
phology of Si/Si0 2 interfaces [39j \M SU [42], making 
it a natural guess for the ID/ED morphologies in our 
study. For a given set of disorder parameters {D\,D2}, 
we compute the ensemble average of N samples to ob- 
tain the conductance a pn , where N is chosen to be from 
100 — 200 samples. The root-mean-square (RMS) and 
auto-correlation length (AL) of the interface roughness 
morphology are also computed. We defined AL to be 
the length at which the cross-correlation is 50% of the 
auto-correlation. It can be shown that AL depends on 
D2 and is relatively insensitive to Di. A larger D2 will 
yield more higher frequency components in I(y), thereby 
decreasing the AL. 

In the presence of disorder, the device Hamiltonian's 
on-site and coupling energies have to be modified accord- 
ingly. For Vi, we have, 



x' < -d 
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x' > d 



(9) 



where x' = x + I pn (y) and epn = \ e F — e° being the 
Dirac point energy. ep P is defined similarly. We had as- 
sumed that the pn junction is linearly graded, where the 
spatial extent of n-p transition (known as the "depletion 
width" ) is denoted by d. I pn is a one dimensional rough- 
ness profile generated using the above procedure. For Uj 
(in eV), we have, 



\t(x,y)\ 



3 , I b {x) < y < W + I t {x) 
, otherwise 



(10) 



I b i t is a one dimensional roughness profile describing the 
line edge roughness for the bottom and top edges respec- 
tively. 

As an illustration, we modeled the magneto-transport 
across an armchair graphene pn junction with/ without 
ID in a magnetic field of B — 10T with energies ep-n = 
O.leV and 6f p = O.lbeV. Fig. [2] plots the longitudi- 
nal current component j x (x : y) of this device for different 
device width. The so-called "snake states" (i.e. current 
density oscillating back and forth the n/p regions) propa- 
gating along the pn interface can be observed [23] [43] [44] , 
which remains prominent even in the presence of ID. The 
snake states terminate when the pn interface meets the 
top edge, where a choice between the paths leading to 
the n or p medium must be made. The ribbon's width, 
which also determines the valley isospins of the first LL 
along each edge, play a pertinent role in deciding which 
path is taken, as will be elaborated upon in the following 
section. 



III. VALLEY ISOSPINS ON THE BALLISTIC 
QUANTUM HALL PLATEAUS 

In this section, we examine the conductance across a pn 
junction when there is no disorder i.e. the ballistic limit. 
The ribbon's chirality plays an important role in deter- 
mining the junction's conductance. Fig. |TJl depicts the 
notation on ribbon's chirality used in this work. In this 
limit, it was shown that conductance plateaus a pn could 
emerged due to the valley isospins [9] and wavefunction 
parity [10] of the ground state Landau level. The former 
effect concerns the armchair edge ribbon, while the latter 
for zigzag edge. An excellent review on this subject has 
been written by Beenakker [23], where we will highlight 
and expand on some of the key results in the remainder 
of this section. 



Valley isospins along the edges: 
Definitions and conventions 



We write the Dirac equation for graphene as, 



Vfp-a 
Vfp-a 



(11) 



where ^ = {^a^b^—^b^a) and VKVO f° r the K(K') 
valley wavefunction. We are interested in ^ along the 
ribon's edges. It is a convenient convention to write ^ 
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FIG. 3: (a) Junction conductance as a function of depletion width 
in the clean limit i.e. no disorder, for armchair ribbons of different 
widths for the case of filling factor (2, 2). (b) Same as (a) excepts 
for filling factor of (6, 2). 



along the edges in the following form, 
9 = (v • f ) <g> (ft • S) 9 



(12) 



where v is the edge valley isospin (for ^ along the 
edges) and r is just the Pauli matrices for the isospin 
part, ft depends on the edge type i.e. ft = (0, 0, 1) for 
zigzag and ft = (±1,0,0) for bottom/top edges of arm- 
chair ribbons [23 j. One can show that Eq. [l2| effectively 
expressed the boundary conditions of the edges . For an 
armchair ribbon, it can be shown that the isospin along 
the top/bottom edges (vt and vb) obeys the following, 



v T • fJ B = cos (AW + 7r) = cos 6 



(13) 



where A = 47r/3a and a is the lattice constant of 
graphene. From Fig. [TJ we have W = a(l + where I 
is the number of carbon layers. See Appendix A for the 
detail algebra. 

Next, one makes the assumption that the ground state 
LL's wavefunction, denoted by |0), could be approxi- 
mated by the edge wavefunctions [9]. This allows one 
to write the wavefunction overlap between the ground 
state LL's wavefunction along the top/bottom edges, 

(Or | B ) « (* T | 9 B ) (14) 
= ^|a| 2 + |fr| 2 ^ (1 + vt ' v B + vt x vb) 

where we denote (^a^b) = (a, b) and (i^a^b) is ob- 
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tained through Eq. 
\b\ 2 = |, we finally arrive at, 



By using the fact that 



B. Valley isospin in armchair ribbon 

Fig. [3] shows the ballistic conductance of armchair 
edge type ribbons as a function of depletion width. The 
n/p regions are biased at e Fn / Fp respectively, and the 
built-in potential (assumed to be linearly graded across 
the junction) is given by e Fn + epp- Fig. [3^l plots the 
conductance for biasing conditions corresponding to the 
Landau filling combinations of (v n ,v p ) = (2,2). Rib- 
bons with different number of carbon layers along the 
width are considered, where the inter-layer separation is 
^/3L [L being the carbon-carbon bond length). These 
ribbons exhibit conductance plateaus of \ and 2 at suf- 
ficiently large depletion width of > 25nm, where typical 
length scale of depletion width in experiments employ- 
ing top/bottom gating scheme are usually several times 
larger than 25nm [45]. These plateaus emerge as long 
as the depletion width is sufficiently large, irregardless of 
the filling factor combinations. 

Fig. [3]o plots the case when (v n ,v p ) = (6,2), and we 
had also checked that these ballistic plateaus remain in- 
tact when (v ni v p ) = (6,6). It is observed that increas- 
ing depletion width filters off the higher Landau levels, 
such that only the zeroth mode Landau edge states con- 
duct through the junction. This is reminiscent of the 
more well-known filtering action of off-normal transverse 
modes by a pn junction in the zero magnetic field case 
[45j|46], although the physics in this context is completely 
different. This might find applications in devices that use 
the Landau levels as an information bit [47]. However, 
pn interface disorder would negate such filtering action, 
to be discussed in Sec. HVl 



Tworzydlo and co-workers [9 attributed the origin of 
the ballistic plateaus to the different valley isospins of the 
th LL at the two edges of the ribbon. The valley isospins 
can be determined from their respective boundary con- 
ditions [23] . When the number of layers satisfy the con- 
dition 3M + 1, where M is an integer, it would exhibit 
conductance plateaus of 2. Note that this is also the same 
criterion for obtaining a metallic armchair ribbon. This 
device is illustrated in Fig. [2^i and c for semiconducting 
and metallic type armchair ribbon respectively. We de- 
fine the following scattering states for the th LL; \0 u b) 
and \0 u t) for the incoming and reflected states, where the 
subscript n/p and T/B denotes the electron/hole medi- 
ums and top/bottom edges respectively. We can compute 
the reflection coefficient by following the prescription in 



r = {0 nT \S\0 nB ) 



|(0t|0b)| 2 «|(1 + CO50) 



(15) 



Vs. 



3M+1 
otherwise 



(16) 



Eq. 15 is a simple result [9] that we will employ in the 



remaining of this section to draw some simple conclusions 
about the magnetotransport properties in a graphene pn 
junction. 



The scattering matrix S describing the evolution of the 
incoming scattering state \0 u b) along the pn interface 
can be simply described by a unit matrix with a constant 
phase factor. The conductance plateaus a pn is then given 
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FIG. 4: Intensity plot of the non-equilibrium electron density, 
logio(n), for various ribbons with width « lOOnra. We plotted 
for the case (a) armchair ribbon (no disorder) with 401 carbon lay- 
ers along the width, (b) zigzag ribbon (no disorder) and (c) zigzag 
ribbon (with interface disorder). The magnetic field is assumed to 
be 10T, while the depletion width 25nm. 

by i, 

on . |2 v / 2 3M+1 
a pn = 2(l-\r\ ) = j i otherw . ge (17) 

This argument requires the assumption that the valley 
isospin obeys the orthogonality identity v u t • v p t — v u b • 
v p b = for current conservation to hold [56] . 

C. Valley isospin in zigzag ribbon 

For zigzag ribbons, a similar width dependent effect 
can be observed, except that the conductance plateaus 
are and 2.0 for zigzag and anti-zigzag ribbons respec- 
tively [9 j. These ballistic plateaus cannot be explained 
by the similar valley isospin argument as used for the 
armchair case. As first pointed out by Akhmerov and 
co-workers [10] , the reflected and transmitted edge states 
both reside on a valley different than the incident state 
i.e. v u t • v p t = v u b • v p b = 1. Therefore, current conser- 
vation entails an inherent intervalley scattering process. 
Fig. [4] shows the intensity plot for the non-equilibrium 
electron density. Local peaks in the electron density 
can be observed at the positions where the pn inter- 
face and the ribbon edge meet. These are signatures 
of the intervalley scattering processes that have taken 



place. Heuristically speaking, one can view the propa- 
gating states along the pn interface as similar to that of 
an armchair edge, where the valley isospin is an equal 
weight superposition of the two valleys. The two times 
scattering process take the Landau state from one valley 
to a superposition and then finally to the other valley. 
In other words, the valley isospin information of the in- 
cident Landau edge state is intrinsically diluted after the 
first scattering process. On a related note, the local den- 
sity of states in quantum Hall regime had been recently 
probed through scanning tunneling microscopy measure- 
ment [48 . It would be of fundamental importance to 
experimentally verify the existence of these localized in- 
tervalley scattering centers. 

Akhmerov and co-workers [10 worked out the trans- 
mission/reflection coefficient for the edge states in the 
ballistic case (no magnetic field case) to be 

r = (0 nT \S\0 nB ) = l I (18) 
N 1 [0 anti-zigzag v J 

where S should embody the valley scattering processes. 
The conductance a pn can then be computed in the same 
fashion as the armchair case. 

f0 zigzag 
apn ~ \ 2 anti-zigzag [ * } 

By defining (a pn ) to be the conductance plateaus av- 
eraged over the ribbon width, we obtained (cr P n) ac — 

(°pn) zz = I- 

IV. TRANSITION TO OHMIC TYPE 
QUANTUM HALL PLATEAUS 

In this section, we examine the role of ID and ED on 
the quantum Hall plateau transition from ballistic-type 
to Ohmic-type i.e. a pn to a pn . Our numerical results 
show that this transition can be ID/ED-mediated de- 
pending on the filling factor combination (v n ,v p ) and 
the ribbon type, as summarized in Fig. [TJd. The ob- 
jective of this section is to present the argument as to 
why valley isospin dilution in general is necessary for the 
ballistic-Ohmic quantum Hall plateau transition, corrob- 
orated with numerical simulation results. 



A. Armchair edge ribbons 

Fig. [5^i studies the junction conductance in the pres- 
ence of pn interface disorder only, for an armchair rib- 
bon of 401 layers. The ballistic plateau for this device is 
& P n — \ since the number of carbon layers is 7^ 3M + 1. 
From a (z/ n , v p ) -independent a pn in the non-disordered 
limit, the junction conductance begins to adopt different 
(z/ n , ^-dependent values as the ID RMS increases. The 
conductance saturates at large enough disorder strength. 
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FIG. 5: (a) Junction conductance as a function of IR disorder characterized by RMS for an armchair ribbon. All devices with open 
symbols have IR disorder of AL w 7.5nm, whereas those with solid symbols have AL w llnra. No ER disorder and a ribbon width of « 
0.1/ira, i.e. 401 carbon layers, is assumed. Bias condition corresponding to filling factor combinations of (2, 2), (2, 6) and (6, 6) are plotted, 
see right inset for details, (b) Junction conductance as a function of ER disorder characterized by RMS (expressed in terms of number of 
layers) for an armchair ribbon. IR disorder assumed to have AL « 7.5nm and RMS ~ llnm. (c) Same as (a), except for zigzag ribbon 
counterpart. In these plots, each data point is an ensemble average over 100 samples. The magnetic field is assumed to be 10T and the 
depletion width 25nm. The dashed lines are plots of (a pn ) ±var(cr pn ). Thick dashed lines are drawn as guide to the eye so as to indicate 
the onset of conductance plateau. 



However, we observed that only the junction conduc- 
tance for the (6, 6) case approaches the Ohmic values 
of cr pn (6,6) = 3. In particular, the junction conduc- 
tance for (2, 2) is extremely robust against IR disorder. 
The junction conductance for (2,6) was enhanced with 
increasing IR disorder. However, it reaches a conduc- 
tance plateau of only w |, lesser than the Ohmic value 
of (T pn (2,6) = §. In computing the conductance for a 
given RMS, we performed an ensemble averaging over 
100 devices with different roughness configurations. Fig. 
[5jb) plots the junction conductance as a function of edge 
disorder RMS, but with a fixed ID. Evidently, edge dis- 
order with only a RMS of one carbon layer would suf- 
fices in inducing the plateau transition from a pn to a pn . 
These observations suggest the following proposition, "In 
an armchair edge ribbon with filling factor combination 
of (2, v p ), the plateau transition from a pn to a pn is both 
ID/ED-mediated. " 

We consider a Chalker-Coddington [14 type argument 
in support of the above proposition. This model con- 
siders the following facts: (a) in the absence of time re- 
versal symmetry, the electronic states exhibit only uni- 
directional transmission (b) the scattering wave func- 
tion follows approximately the equipotential lines of the 
random potential, which is shown in Fig. [2] Consider 
(vn,v p ) = (2, 6), the scattering state for a particular spin 
along the pn interface can be expressed as, 



four modes scattering can be parameterized as, 



|*i> = c Kb) + ci |0 P > + c 2 |l p ) + c 3 \l p ) 



(20) 



where |0 P ), and \l' p ) are the ground and first ex- 
cited states of the LL in the p medium respectively. We 
have &i = (1, 0, 0, 0) at the beginning of the pn interface. 
We can define a "saddle point" to be where two Landau 
modes i and j undergo mode mixing, characterized by a 
scattering matrix which evolves the scattering state 
in a unitary manner. The effective unitary matrix for 



c 2 


sc 


s 2 


—sc 


—sc 


c 2 


sc 


s 2 


s 2 


—sc 


c 2 


sc 


sc 


s 2 


—sc 


c 2 



(21) 



where s = sin(/3) and c = cos(/3). As usual, the ac- 
companied phase factors is implicit [14]. The parameter 
f3 characterized the degree of mode-mixing i.e. (3 = 0, j 
denotes minimum/maximum mixing. Undergoing a suffi- 
cient amount of mode mixing processes 5, the wavef unc- 
tion at the end of the pn interface could be expressed as, 



|*/) = S(p 1 )S(f3 2 )S(f3 3 )...\^ i ) (22) 
\0ub) + |0 P > + \c^ \\ v ) + \c^ |i;> 



2 e 



The final state is then said to have "completely equi- 
librated". The reflection probability can then be com- 
puted, 



1(0, 



I* 



\\(0uT |O„B>| 2 |e 00 



(23) 



In arriving at the above result, we had make use of the 
orthogonality relation (0 u t = 0. We also assumed 
that 1 P ) retains the isospin information of the incident 
scattering state and therefore yielding us (0 u t| p ) = 
(0tit| n #). By making use of the fact that the phase 
term averaging over a sufficiently large ensemble yields, 



O 0O 



1 



ensemble 



(24) 



the junction conductance (including spin) at filling factor 
(2, 6) could then be expressed as, 

(7p n (2,6) * 2 (l - \ |(0 nT |0n£>| 2 ) (25) 
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FIG. 6: (a) depicts the theoretical Ohmic plateaus cr pn (v n , u p ) [8] as a function of 711/712 where the different colors represent the filling 
factors, (b) — (d) plots the linescan for the following cases (i) v n = v v (ii) v n = 2 and (iii) v n = 6 respectively, and compare the simulation 
results with experimental data. The experimental data are taken from 7 , for a two terminal pn junction quantum Hall measurement at 
B = AT. The simulations are done at B = 10T. In the experimental data, the gate oxide capacitance is used as a fitting parameter. 



For the case where the number of carbon layers of the 
armchair ribbon / 3M + 1, Eq. [25] yields a pn (2,6) = |. 
Repeating the above analysis formling factor (2,2), we 
can show that <r pn (2,2) = <j pn (2,2). The results from a 
simple Chalker-Coddington analysis are in excellent cor- 
roboration with what we obtained numerically from nu- 
merical calculations as shown in Fig. [5^. 

Remarkably, the above analysis predicts that the junc- 
tion conductance for an armchair ribbon with 3M + 1 
carbon layers will always be perfectly conducting i.e. 
<r pn (2,6) = <j pn (2,2) = 2. Fig. [2^-d depicts the spatial 
current density for this device at (2,6) filling factor. In- 
deed, the transmission remains perfect in the presence of 
ID, despite the fact that mode mixing via ID has taken 
place. This fact unequivocally demonstrates that elec- 
tron and hole Landau modes mixing via ID alone is not 
sufficient in achieving the ballistic- Ohmic plateau tran- 
sition, at least for armchair ribbons. It is evident from 
Eq. [25] that the Ohmic result can be obtained if and only 
if |(0 n T \^ub)\ 2 = §• This is only possible if the isospin 
information on the T/B edges is completely diluted, e.g. 
via edge disorder. 



B. Zigzag edge ribbons 

Fig. (5)3 study the junction conductance in the presence 
of pn interface disorder only, for a zigzag ribbon. The bal- 
listic plateau for this device is a pn = 0. In this case, the 
d pn — > (jpn transition in the presence of ID is succinctly 
illustrated. Defining Dq to be the disorder-related length 
scale where the junction conductance begins to plateau 
off, we obtained Do ~ 3,6,8nm for the (2,2), (2,6) and 
(6, 6) case respectively. Similar Do are observed for the 
armchair counterpart devices in Fig. [5^. This leads us 
to the following conclusion, "In a zigzag edge ribbon, ID 
alone is sufficient for the plateau transition a pn to <J vn . 
In addition, it requires about the same disorder strength 
as its armchair counterpart to obtain the Ohmic plateaus, 
albeit ED is required for the latter. " 

Repeating the Chalker-Coddington type argument, the 



reflection probability for the (2, 6) filling factor case can 
be written as, 



(26) 



18 and that (0 nT | P ) = 1 



by making use of Eq. 
(0tit| n #). This then gives us a junction conductance 



of <r pn (2,6) 
at <7p n (2, 2) 



= <Xp n (2,6). Similarly, we can arrive 
<7 pn (2,2). Coincidentally, the signa- 



tures of the parity effect coming from the ground state 
LL is completely neutralized by the simple fact that 
(0 n T\ P ) = 1 — (0 n T\ u b)- A similar physical picture 
applies for the anti-zigzag ribbon case. 



V. BENCHMARKING WITH EXPERIMENTS 

In this section, we benchmark our numerical results 
against the two terminal quantum Hall measurement per- 
formed by Williams et al. [7 j. We should emphasize that 
bulk disorder was not included in present numerical sim- 
ulation, therefore additional longitudinal resistance con- 
tributions which might exist in actual experiments are 
not captured [49 . In the experiment, a top/bottom gate 
is employed to bias the two junctions, allowing control 
over the sign/magnitude of the charge density (711/712) 
residing in each junction. Fig. |6^i depicts the theo- 
retical Ohmic plateaus cr pn (v ni v p ) [8 as a function of 
the electron density n\ / (the different colors represent 
the filling factors). In the numerical calculations, the 
electron density is obtained by taking the trace of the 
electron correlation function G n (e j p), an energy resolved 
quantity. This is defined to be G n = GE^G*" where 
S in = -2Im(E a + E d ) for T = OK. Electron density 
in either the n/p medium can then be computed via the 
integral n = J (G n ) de, where the averaging (...) is per- 
formed over the spatial dimension. 

Fig. j6|>d plots the linescan for the following cases (i) 
v n — v p (ii) v n = 2 and (iii) v n = 6 respectively. In gen- 
eral, the numerical results show satisfactory agreement 
with the experiments. As previously addressed [7 , the 
junction conductance with lower filling factors such as 
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cr pn (2, 2) and <r pn (2, 6) plateau off at the expected Ohmic 
values of 1 and | respectively. However, higher plateaus 
such as (6,6), (6,10) could not be observed experimen- 
tally. This suggests that the interface disorder in the 
experiment is smaller than that necessary for complete 
Landau mode mixing of the higher plateaus. Fig. [5^i 
indicates that ID with RMS and AL of « lOnm should 
be sufficient. A more recent experiment by Lohmann et 
al. [50] which employed chemical doping methods to cre- 
ate pn junctions should exhibits a larger ID. Although 
higher plateau measurements were not done in their ex- 
periments, their lowest plateau <r pn (2,2) exhibits a more 
precise plateau than that reported in [7j. However, the 
improved precision in measurement is also a direct result 
of four terminal measurement [50 . 

The general characteristics of the conductance dur- 
ing plateau transition also agree qualitatively with that 
of experiments i.e. the decreasing slope da pn /dn for 
the higher plateaus transitions. This is attributed to 
the smaller peak in the density of states (of the so- 
called extended states) of the higher LL [51]. Another 
general remark can be made about the conductance 
plateau. Theoretically, the inter-LL energy spacing de- 
creases with higher LL in a manner that is proportional 
to y/n — \Jn — 1, n being the LL's index. However, the 
inter-LL spacing as function of charge density n x / 2 (inter- 
LL spacing herein denoted as S n ) is approximately equi- 
distant as shown in Fig. [6^. The conductance plateau 
width for the filling factor v n = v p could then be ex- 
pressed as 5 n — r]/(d<j pn /dn), where r] = 2e 2 /h. Since 
d(jp n /dn decreases with increasing the plateau width 
has to also decrease accordingly. This renders the obser- 
vation of conductance plateau at higher LL more chal- 
lenging in experiments. Nevertheless, our study suggests 
that the observations of junction plateaus up to filling 
factor of 6 might be possible. 



VI. SUMMARY 

In summary, we conducted a systematic study of 
graphene pn junction conductance in the quantum Hall 
regime. Often, the disorder in magneto-transport cal- 
culations are modeled in an implicit manner through 
a random on-site energy fluctuation. In this work, we 
had undertaken the effort to explicitly modeled the vari- 
ous disorder in order to uncover the underlying physical 
mechanisms played by pn interface and edge disorders on 
the ballistic-Ohmic quantum Hall plateau transition. We 
found that the former mechanism equilibrates the elec- 
tron/hole Landau modes along the interface while the 
latter dilute the isospin information of the ground state 
LL in an armchair edge ribbon. For a given Landau filling 
factor combination, (y n ,v v \ we found that both zigzag 
and armchair ribbons require about the same pn inter- 
face disorder strength to recover the Ohmic plateaus, al- 
beit edge disorder is required for the latter. From our 
numerical calculations, we found that pn interface disor- 
der with a root-mean-square roughness of lOnm is suffi- 
cient in achieving complete mixing of the electron/hole 



Landau modes. However, the sloppiness in the quantum 
Hall plateau transition induced by the pn interface disor- 
der and the accompanied decreasing plateau widths with 
increasing filling factor makes it challenging to observe 
these higher plateaus experimentally. 

From a theoretical standpoint, we argued that the 
mixing of electron/hole Landau modes along the inter- 
face alone does not guarantee the recovery of the Ohmic 
plateaus. We extended the valley isospin argument pro- 
posed by Tworzydlo and co-workers [9] with a Chalker- 
Coddington [14] type argument to highlight this point, 
and corroborated the calculations with numerical simu- 
lations. It is demonstrated numerically that both inter- 
face and edge roughness (or intervalley scattering) are in 
general necessary for the crossover between the two theo- 
retical limits i.e. ballistic and Ohmic. Last, but not least, 
this work underscores the importance of including both 
interface and edge disorder in the modeling of quantum 
Hall transport phenomena, especially when interpreting 
experimental data. 
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APPENDIX A: BOUNDARY CONDITIONS AND 
EDGE ISOSPINS FOR ARMCHAIR RIBBONS 

We consider an armchair ribbon where the two edges 
are at y = ut (top) and y = ys (bottom). Along y = y?, 
the wavefunction \I/ = (ipA,ipB, —^b^a) must satisfy 
the boundary conditions [52] , 



(Al) 
(A2) 



where A = 47r/3a and a is the lattice constant of 
graphene. We can rewrite Eq. A2 in the form ^ = M\I/, 
where, 



(A3) 

(A4) 
and 



- e -i^VT - 




"o r 


e iAyT 


<8> 


1 



M = 

With some matrix algebra we can show that, 
M — (v T • t) (Ht • 5) 



by defining VT=(cos(AyT), 5w(A^),0) 
nT=(— 1, 0, 0). Repeating this procedure for ys, we re- 
quire VB=(—cos(AyB),-sin(AyB),0) and nB = (l,0,0). 
It is straightforward to see that, 



v T • v B = cos (AW + 7r) = cos(0) 
where W = yr — Vb- 



(A5) 



10 



[1] K.S.Novoselov, A.K.Geim, S.V.Morozov, D.Jiang, 

Y.Zhang, S.V.Dubonos, I.V.Grigorieva, and A.A.Firsov, 

Science 306, 666 (2004). 
[2] G.W.Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 
[3] K.S.Novoselov, A.K.Geim, S.V.Morozov, D.Jiang, 

M.I.Katsnelson, I.V.Grigorieva, S.V.Dubonos, and 

A.A.Firsov, Nature 438, 197 (2005). 
[4] Y.Zhang, Y.W.Tan, and P.Kim, Nature 438, 201 (2005). 
[5] C.W.J.Beenakker, C.Emary, M.Kindermann, and 

J.L.v.Velsen, Phys. Rev. Lett. 91, 147901 (2003). 
[6] D.A.Abanin, P.A.Lee, and L.S.Levitov, Phys. Rev. Lett. 

96, 176803 (2006). 
[7] J.R.Williams, L.DiCarlo, and C.M.Marcus, Science 317, 

638 (2007). 

[8] D.A.Abanin and L.S.Levitov, Science 317, 641 (2007). 
[9] J.Tworzydlo, I.Snyman, A.R.Akhmerov, and 

C.W.J.Beenakker, Phys. Rev. B 76, 035411 (2007). 
[10] A.R.Akhmerov, J.H.Bardarson, A.Rycerz, and 

C.W.J.Beenakker, Phys. Rev. B 77, 205416 (2008). 
[11] W.Long, Q.F.Sun, and J.Wang, Phys. Rev. Lett. 101, 

166806 (2008). 
[12] J.Li and S.Q.Shen, Phys. Rev. B 78, 205308 (2008). 
[13] P.W.Anderson, Phys. Rev. 109, 1492 (1958). 
[14] J.T.Chalker and P.D.Coddington, J. Phys. C 21, 2665 

(1988). 

[15] B.Huckestein, Rev. Mod. Phys. 67, 357 (1995). 

[16] W.R.Frensley, Rev. Mod. Phys. 62, 745 (1990). 

[17] R.Landauer, Philos. Mag. 21, 863 (1970). 

[18] M.Buttiker, Phys. Rev. Lett. 57, 1761 (1986). 

[19] S.Datta, Electronic Transport in Mesoscopic System, 

Cambridge University Press (1997). 
[20] H.Haug and A.P.Jauho, Springer Series in Solid State 

Sciences (Springer, New York) 123 (1996). 
[21] C.S.Lent and D.J.Kirkner, J. Appl. Phys. 67, 6353 

(1990). 

[22] A. Neto, F.Guinea, N.M.R.Peres, K.S.Novoselov, and 
A.K.Geim, Rev. Mod. Phys. 81, 109 (2009). 

[23] C.W.J.Beenakker, Rev. Mod. Phys. 80, 1337 (2008). 

[24] P.R.Wallace, Phys. Rev. 71, 622 (1947). 

[25] R.Saito, G.Dresselhaus, and M.S.Dresselhaus, Physi- 
cal Properties of Carbon Nanotubes (Imperial, London) 
(1998). 

[26] K.Wakabayashi, M.Fujita, H.Ajiki, and M.Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

[27] Y.Imry, Introduction to Mesoscopic Physics (Oxford Uni- 
versity Press, New York) (2002). 

[28] G.D.Mahan, Many Particle Physics (Plenum, New York) 
(1990). 

[29] M.P.L.Sancho and J.M.L.Sancho, J. Phys. F: Met. Phys. 

14, 1205 (1984). 
[30] S.Nonoyama and A.Oguri, Phys. Rev. B 57, 8797 (1998). 
[31] M.P.Anantram, M.S.Lundstrom, and D.E.Nikonov, 

Proc. IEEE 96, 1511 (2008). 
[32] G.Grosso, S.Moroni, and G.P.Parravicini, Phys. Rev. B 

40, 328 (1989). 
[33] T.Low and J.Appenzeller, Phys. Rev. B (2009). 



[34] D. Fisher and P.A.Lee, Phys. Rev. B 23, 6851 (1981). 
[35] E.H.Hwang, S.Adam, and S.D.Sarma, Phys. Rev. Lett. 

98, 186806 (2007). 
[36] I.F.Herbut, V.Juricic, and O.Vafek, Phys. Rev. Lett. 

100, 046403 (2008). 
[37] J.Martin, N.Akerman, G.Ulbricht, T.Lohmann, 

J.H.Smet, K.V.Klitzing, and A.Yacoby, Nat. Phys. 

4,144 (2008). 

[38] A.K.Gupta, T.J.Russin, H.R.Gutie, and P.C.Eklund, 

Nano Lett. 3, 45 (2009). 
[39] S.M.Goodnick, D.K.Ferry, C.W.Wilmsen, Z.Liliental, 

D.Fathy, and O.L.Krivanek, Phys. Rev. B 32, 8171 

(1985). 

[40] T.Low, M.F.Li, W.J.Fan, S.T.Ng, Y.C.Yeo, C.Zhu, 

A. Chin, L.Chan, and D.L.Kwong, IEEE Int. Elec. Dev. 

Meet. p. 151 (2004). 
[41] T.Low, M.F.Li, G.Samudra, Y.C.Yeo, C.X.Zhu, A.Chin, 

and D.L.Kwong, Trans. Elec. Dev. 52, 2430 (2005). 
[42] T.Low and D.Ansari, Phys. Rev. B 78, 165301 (2008). 
[43] J.M.Pereira, F.M.Peeters, and P.Vasilopoulos, Phys. 

Rev. B 75, 125433 (2007). 
[44] V.Lukose, R.Shankar, and G.Baskaran, Phys. Rev. Lett. 

98, 116802 (2007). 
[45] T.Low, S.Hong, J.Appenzeller, S.Datta, and 

M.S.Lundstrom, Trans. Elec. Dev. 56, 1292 (2009). 
[46] V.V.Cheianov and V.I.Fal'ko, Phys. Rev. B 74, 041403 

(2006). 

[47] C.W.J.Beenakker, Proc. Int. School Phys. E. Fermi 162 
(2006). 

[48] K.Hashimoto, C.Sohrmann, J.Wiebe, T.Inaoka, 
F.Meier, Y.Hirayama, R.A.Romer, R.Wiesendanger, and 
M.Morgenstern, Phys. Rev. Lett. 101, 256802 (2008). 

[49] J.R.Williams, D.A.Abanin, L. DiCarlo, L.S.Levitov, and 
C.M.Marcus, Phys. Rev. B 80, 045408 (2009). 

[50] T.Lohmann, K.V.Klitzing, and J.H.Smet, 

http: / /arxiv.org/abs /0903.5430 (2009) . 

[51] P.Streda, J. Phys. C 15, L717 (1982). 

[52] L.Brey and H.A.Fertig, Phys. Rev. B 73, 235411 (2006). 

[53] K.H.Bevan, T.Low, and H.Guo, J. Appl. Phys. 105, 
093709 (2009). 

[54] A more recent experiment by Lohmann et al. [50], which 
employed a four terminal measurement, observed promi- 
nent plateaus for only (2,2). The higher filling factor 
combination values were not reported. In their experi- 
ments, the pn junctions were created via chemical dop- 
ing, instead of the usual top/bottom gating scheme [7]. 

[55] However, one could easily extends the simple tight bind- 
ing Hamiltonian to include higher orbitals by extracting 
the orbital coupling energies from a density functional 
calculation as described in [53] . 

[56] The sum of transmission and reflection probability 
must be unity. In other words, we should have 1 — 
\(0 nB \ nT )\ 2 = \(0 nB \ p t)\ 2 . This only holds true if 

VnT • VpT — VnB ' V p B = 0. 



